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The coupling of antimatter to gravity is of general interest because of conceivable cosmological 
consequences ("surprises") related to dark energy and the cosmological constant. Here, we revisit the 
derivation of the gravitationally coupled Dirac equation and find that the prefactor of a result given 
previously in [D. R. Brill and J. A. Wheeler, Rev. Mod. Phys. 29, 465 (1957)] for the afiine connection 
matrix is in need of a correction. We also discuss the conversion the curved-space Dirac equation from 
East-Coast to West-Coast conventions, in order to bring the gravitationally coupled Dirac equation 
to a form where it can easily be unified with the electromagnetic coupling as it is commonly used in 
modern particle physics calculations. The Dirac equation describes anti-particles as negative-energy 
states. We find a symmetry of the gravitationally coupled Dirac equation, which connects particle 
and antiparticle solutions for a general space-time metric of the Schwarzschild type and implies 
that particles and antiparticles experience the same coupling to the gravitational field, including 
all relativistic quantum corrections of motion. Our results demonstrate the consistency of quantum 
mechanics with general relativity and imply that a conceivable difference of gravitational interaction 
of hydrogen and antihydrogen should directly be attributed to a a "fifth force" ("quintessence"). 
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I. INTRODUCTION 

In view of the recent dramatic progress of antimatter 
gravity experiments [1, 2], it seems indicated to reex- 
amine the theoretical status of antimatter coupling to 
gravity. A number of experimental collaborations are ac- 
tively pursuing related experiments [3-8] . A key factor in 
recent experimental progress [2] of the ALPHA collabo- 
ration has been their special Penning-Ioffe trap which 
simultaneously traps both positrons as well as antipro- 
tons. Superimposed on the Penning trap fields (which 
trap the charged constituent particles), the antihydrogen 
atom Ioffe trap employed by ALPHA relies on a strong 
octupolc magnetic field configuration generated by eight 
superconducting current bars, which wind back on them- 
selves in a snakc-likc pattern, glued to the inner chamber 
of the ALPHA experiment by a three-dimensional wind- 
ing machine at Brookhaven National Laboratory (see 
Ref . [9] ) . This leads to an effective trapping of positrons 
and antiprotons, and antihydrogen atoms. 

Antimatter gravity experiments aim to test the inter- 
action of antihydrogen atoms with gravitational fields. 
According to general relativity [10], gravitational inter- 
actions can be described by the induced space-time cur- 
vature around massive objects. Furthermore, on the clas- 
sical level, the motion of a particle in curved space-time 
is described by the following geodesic equation [10] 

d 2 X P „ dxP dx a 

which implies that a particle of mass m experiences 
a "force" = md 2 x fl /ds 2 and moves along a zero 
geodesic in the gravitationally curved space-time (s is 
the proper time). Here, the T* 1 pa are the Christoffel sym- 
bols [10], derived from the curved-space metric as 



follows, 

apa ~ 2 V dxP + dx° dx a J ' [ ' 

where the Einstein summation convention is used. The 
T 11 pa arc derived from the Y apa by raising the first in- 
dex with the help of the metric, i.e., F M pa = r ap „. 
We should clarify that in the current article, in a some- 
what non-standard notation, the symbol <?„„ will be re- 
served for the flat-space metric in the following, whereas 
~g >J,a denotes the metric of curved space. If, according to 
Einstein's equivalence principle, we assume that gravita- 
tional mass and inertial mass are proportional to each 
other, then classical gcomctrodynamics [10], on the basis 
of Eq. (2), makes the unique prediction that the force on 
a particle and antiparticle in a gravitational field are the 
same, provided the mass of particle and antiparticle are 
equal, i.e., both particle as well as antiparticle motion are 
described by Eq. (1). However, on the quantum level, the 
situation is less clear. It is often argued [7] that "general 
relativity is incompatible with quantum mechanics" and 
that, assuming rather peculiar couplings of antimatter 
to gravity [11], one can imagine that antimatter actually 
is repulsed by gravity. This observation provides part of 
the motivation for a number of antimatter gravity exper- 
iments currently under preparation [3-8]. 

Here, we reexamine the status of theoretical predic- 
tions regarding the coupling of Dirac particles and an- 
tiparticles to curved space-time. Indeed, closer inspec- 
tion shows that considerable insight into the gravitational 
coupling of antiparticles can be gained based on rather 
straightforward generalizations of previous treatments 
which rely on a combination of relativistic quantum me- 
chanics with general relativity. We note the works of Brill 
and Wheeler [12], Boulwarc [13], and Soffcl, Mullcr and 
Greincr [14]. The Dirac equation [15, 16] describes both 



2 



particles and antiparticlcs simultaneously, and symme- 
tries of the solutions which connect particles and antipar- 
ticles are therefore relevant for antigravity experiments. 
The equations derived from the formalism of Refs. [12- 
14] have never been analyzed in terms of their conse- 
quences for antimatter coupling to gravity, to the best 
of our knowledge. We find that it is highly indicated to 
revisit a number of aspects of the derivation. We employ 
units with H = c = eg = 1. 



II. FORMALISM 

Antihydrogen consists of two spin-1/2 particles, the 
electron and the proton. Spin-1/2 particles are described 
by the Dirac equation. In curved and flat space-time, 
respectively, the anticommutators {•,•} of the Dirac 7 
matrices fulfill the algebraic relations 



(3) 



where the curved-space metric is ~g ilv with fi, v = 0,1,2,3, 
while the local flat-space metric ( "vierbein" ) in our con- 
ventions is ~g^ v = diag(l, — 1, — 1, —1). The precise form 
of the 7 M (x) matrices depends on the space-time geom- 
etry and in particular, on the space-time coordinate x. 
We use the "West-Coast" signature diag(l, — 1, — 1, — 1) 
for the free-space metric g instead of the "East-Coast" 
conventions diag(l, 1, 1, —1) or diag(— 1, 1, 1, 1), in order 
to ensure compatibility with the sign convention usually 
adopted in the modern particle physics literature [17-20]. 

This is different from the sign conventions adopted in 
the traditional literature on general relativity [10], and 
also different from the sign conventions used in previ- 
ous works on the gravitationally coupled Dirac equa- 
tion [12-14]. E.g., the conventions of Misner, Thorne 
and Wheeler [10] are given in Eq. (2.10) on page 53 
of Ref. [10] and involve a metric diag(— 1, 1, 1, 1). with 
/i, v = 0,1,2,3. If we ever wish to study the combined 
"gravito-magnetic" effect of gravitational and electro- 
magnetic fields on Dirac particles simultaneously, and 
conceivably use established results for the electromag- 
netic sector in a perturbative, then it is helpful to convert 
the gravitational Dirac equation into "West-Coast" con- 
ventions, because these are used in the particle physics 
and quantum electrodynamics (QED) literature. 

One might ask why we are using the tilde in order 
to denote the flat-space metric, not just the plain sym- 
bol g^. The answer to that question is as follows. We 
would like to be as unique in our notation as possible, and 
avoid possible confusion upon comparison with the liter- 
ature [12-14]. In Refs. [12-14], the curved-space Dirac 
matrices are denoted as 7 M , but in the particle physics 
literature [17-20], one denotes the flat-space matrices as 
7**. There is no way to unify the notations without intro- 
ducing some ambiguity, and we have therefore decided 
to differentiate the matrices either by overlining or using 
the tilde, making their identification unique. 



The flat-space action for the free Dirac particle in spe- 
cial relativity reads as 



d xip(x) (\^d^ — m ) ip{. x ) 



d 4 x^(x) [if^-m] i>{x), 



where ip(x) = ip(x) + a is the Dirac adjoint, = d/dx^ 
is the derivative with respect to x^, and the symmetric 
derivative operator acts as 



±7$ 



A(x) d „ B(x) = A(x) B(x) - B{x) A(x) . (5) 
Furthermore, a is a Hcrmitizing matrix with the property 
a a = T- (6) 



Here, b + denotes the Hermitian conjugate of a matrix b. 
An infinitesimal global Lorentz transformation A and the 
corresponding spinor Lorentz transformation 5(A) in flat 
space then read as 



(7a) 

W , (7b) 
(7c) 

-a/3 
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Here, u/ 



• u;„ — 0. In formulating the generators a 
of spinor Lorentz transformations, we follow the conven- 
tions of Chap. 2 of Ref. [19]. Furthermore, in view of the 
relation 



= 2ig^ 7 — 2 1 cf 7 



(8) 



the 7 matrices are shape-invariant under Lorentz trans- 
formations, 



7'" = A"„ 5(A) 7" 5(A)- 1 =7^ 



(9) 



and ip transforms with the inverse Lorentz transforma- 
tion, ip (x') = ip(x) 5(A)" 1 . This can be shown easily by 
observing that 



a {S(A)) + a = S(A)- 1 . 



(10) 



Standard representations of the flat-space Dirac matrices 
7 include the Dirac and the Majorana representation [19, 
20]. 

The generalization of the Dirac action (4) to curved 
space-time involves two steps: (i) an obvious generaliza- 
tion of the anticommutator relations (3) to curved space, 
{7 M (x),7 M (a;)} = 2g flIJ (x), and (ii) a coupling of the 
derivative operator in the Dirac equation to the gravi- 
tational field on the basis of a covariant derivative, in the 
sense of the replacement 9„ — > V« = — , where V M 
is the covariant derivative and T M = T^(x) is the affinc 
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connection matrix. The action for the Dirac particle in 
curved space-time then reads as 



S = / d 4 x V- det g ip(x) {i^(x)Vp - m) ip(x) , 



(11) 

where det g = det g^ v < is the determinant of the 

space-time metric, and tp = ip + a(x) is the curved-space 
Dirac adjoint, where a{x) is a Hermitizing matrix with 
the local property 

o(ar) a(x) = *f(x) , (12a) 

a(x) (S(L(x))) + a{x) = S{L{x))- 1 . (12b) 

Here, S(L(x)) is the spinor transformation corresponding 
to an infinitesimal, local Lorentz transformation L(x) and 
reads as 



(13a) 



S(L(x)) = l--* af> (x)W afl (x), (13b) 

* ap {x) = \ [rw.fw] • (i3c) 

These equations generalize Eq. (7) to curved space-time 
and ensure that the 7 matrices are shape-invariant under 
Lorentz transformations, 

^(x)=L(xr v S(L(x))r(x)S(L(x))- 1 =r(x). 

(14) 

From now on, we shall suppress the space-time coordinate 
argument x in the 7 and a matrices. The generalization 
of Eq. (8) to general relativity is given by 



[7", a a0 ] = 2 ig" a 7^ - 2 ig^ 7° 



(15) 



One can show this relation using Eq. (3) only. For tensors, 
the covariant derivative V p is well established [see Ex- 
ercise 8.4 on page 211 of Rcf. [10]], but for spinors, a 
nontrivial generalization is required. Let us assume the 
structure 



(16) 



where the affine connection matrix T u remains to be 
determined. We postulate that the covariant derivative 
operator commutes with the current matrix 7 M (x), i.e., 
[7 M (a:), V,] -0. Then, 

7» V^(x) = V„ (l^x)4>(x)) , (17) 
and we can symmetrize Eq. (11) as follows, 

S = J d 4 x v^dct^ ^(ar) \JfP(x) V p - m\ ^{x) . 

(18) 

By variation, the gravitationally coupled Dirac equation 
is obtained as 



(i^V^ -rn)ip(x) =0. 



(19) 



An additional electromagnetic field could be incorporated 
by the replacement V M — > V M + i q A p , where is the 
vector potential and q the charge. However, the gravi- 
tational Dirac equation is primarily interesting when all 
electromagnetic are compensated and the residual grav- 
itational interaction dominates the kinematics. 

The affine connection matrices remain to be deter- 
mined. Using the ansatz (16), one can write the condition 
given in Eq. (17) as follows [12-14, 21, 22], 



0. 



(20) 



The r p /J!y 7 p are the Christoffcl symbols defined in Eq. (2). 
For a Lorentz vector T M , we recall that [10] 



^ — d p T a r a p T\ 
V M T Q = dpT a + r% A T A 



(21a) 
(21b) 



The third and forth term on the right-hand side in 
Eq. (20) represent the spinor structure contributions to 
the covariant derivative of the 7,, matrix. 

' P 

The condition (20) defines the T u matrix up to a mul- 
tiple of the unit matrix. In the vierbein formalism, we 
can represent the 7" matrices in terms of the vierbein 7^ 
matrices as follows, 



Ip = bp" 7a , 



lp — u P la 1 7p a P la 1 

— a a — p — a j a — p 

7 = a P 7 , 7 = bp Y ■ 



(22a) 
(22b) 



pa 



The metric is recovered as 

{%, 7.} = b p a b/ = 2g ap b p a b/=2] 

(23a) 

= a" a a°p {T, 1^} = 2fl a " a°p = ■ 

(23b) 

We note that the matrix with components ~g ( " T is the in- 
verse of the matrix with components g prT , where the en- 
tries of ]j a p and ~g a ^ are identical. It is possible to show, 
using rather lengthy algebra, that the following affine 
connection matrix, 



" 4 9 pa 



db v 







dxP 



ra 
1 



(24) 



with = i [f , 7"], fulfills Eq. (20) for a general metric 
~gP v . Our result (24) differs from the result given in Eq. (8) 
of Ref. [12] in the correction of an obvious and in some 
sense trivial typographical error. Namely, the expression 
T a u i in Eq. (8) of Ref. [12] should be replaced by the 
expression r Q „k, (we have used p for the corresponding 
index of T p , not Tk as in Rcf. [12]). It is less trivial to see 
that a prefactor 1/4 is missing from Eq. (8) of Ref. [12] 
and needs to be supplemented as given in Eq. (24). In 
order to clarify the matter, we should also point out that 
the additional imaginary unit in the prefactor is entirely 
due to our different conventions for the 7 matrices and 
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the flat-space metric which follow modern "West-Coast" 
conventions [19, 20]. 

It is rather lengthy but straightforward to show that 
Eq. (24) solves Eq. (20). One needs to use Eqs. (22) 
and (15) repeatedly, and one needs to observe that the b 
matrix is the inverse of the a matrix, i.e., a k a b p a = 5 p , 
where S is the Kronecker symbol. Furthermore, the re- 
lation T^ ap + = gP a d p g aa is useful in intermedi- 
ate steps of the calculation. Here, p = g^ a T aap with 
r CT Q P given in Eq. (2). 

Using the identity ~cf^ v = i~g^ v — vj V7 y fJ ', it is possible 
to rewrite Eq. (24) in a simpler form, 



V 
4 



l r 



■ (5 P 7„-r%7„) +A P 1 M , (25a) 
(d p bj) a a p - (d p g a0 ) g a 



(25b) 



For a diagonal structure of the metric tensor (the only 
nonvanishing elements are the g a p with a = j3), with 

bj = y[\3^\ and a a p = = l/y^~~~j, the addi- 

tional term A p vanishes. This is the case for the (gener- 
alized) Schwarzschild metric to be discussed below. We 
have checked that, up to the term A p and up to the 
conversion from East-Coast to West-Coast conventions, 
the result (25a) is formally identical to the result pre- 
viously given in Eq. (9) of Ref. [14]. As a byproduct of 
our calculation of the A p , we thus show that the two 
different results for the affine connection matrix given 
in Refs. [12, 14] are equivalent up to the term A p l4 X 4, 
which is proportional to the unit matrix and not deter- 
mined by the defining Eq. (20). For a diagonal metric, 
A p vanishes. To the best of our knowledge, the precise 
form of the term A p has not yet been indicated in the 
literature. 

Our construction of the spinor Lorentz transformation 
in curved space [see Eq. (13)] follows ideas outlined in 
Ref. [13]. However, our result for the covariant deriva- 
tive of a spinor manifestly contains additional terms as 
compared to the result given in Eq. (2.8) of Ref. [13]. 
In particular, in view of the condition (20), it is clear 
that the derivative terms [the first terms in round brack- 
ets on the right-hand sides of Eqs. (24) and (25a)] are 
an essential contribution to the affinc connection matrix; 
these terms seem to be missing in the vierbcin formal- 
ism formulated in later steps of the derivation leading to 
Eq. (2.8) of Ref. [13]. 



III. SCHWARZSCHILD-TYPE METRIC 

A. Radially Dependent Metrics 

In the following, we shall describe an important appli- 
cation of the above formalism which has not yet been 
described in the literature, to the best of our knowl- 
edge. Namely, we shall discuss a generalization of the 



Schwarzschild metric, which describes, to good approxi- 
mation, the gravitational field of a planet, e.g., the Earth. 
In "West-Coast" conventions (where the sign of the time- 
like component is positive), the Schwarzschild metric 
reads as follows [12-14, 23], 



g = diag (e y ', -c A , -r 2 , -r 2 sin 2 6 

j- I 2 1 2 2-2 

= diag u , r,-r ,—r sin 

\ u 

U W=e- = l-^£ 



r s = 2Mq = 2GM f 



(26) 

(27) 
(28) 



Here, the Schwarzschild radius is r s = 2Mq = 
2GM P /c 2 w 0.0089 m, and G is Newton's gravitational 
constant (here, we supplement the factor c~ 2 although 
we use units with c = 1 in this article otherwise). Fur- 
thermore, Mp is the mass of the Earth (or, of the planet 
under consideration). The invariant line element ds 2 in 
the Schwarzschild geometry is given by 



ds 2 = 



1 _ LI dt 2 - 1 - ^ 
r ' V r 

2 I Aa2 i „:„6» j,„2\ 



dr 2 



sin" d(p z 



(29) 
(30) 



The Schwarzschild metric is valid for a spherically sym- 
metric geometry of space. However, it has a problem. 
Namely, as pointed out by Eddington [24] , for the original 
Schwarzschild metric, the speed of light in the radial di- 
rection is not equal to the speed of light in the transverse 
directions; the prefactor in front of the "angular" term 
r 2 (d(9 2 + sin 2 dip 2 ) is not the same as the one in front 
of the "radial term" proportional to dr 2 . This structure 
implies that one has to resort to a highly nonstandard 
representation of the Dirac algebra [14] if one would like 
to separate the gravitationally coupled Dirac equation in 
the original form of the Schwarzschild metric (26). 

For example, without explicit mention, a representa- 
tion of the following form has apparently been used in 
Ref. [14], 



/ -11-2x2 







-1 2 X2 I ' 7 V -112X2 



-il 2 



~2 








' 




7 = 




j ' 




-a 2 


1) 



~5 • ~0 ~\ ~2 ~3 

7 =17777 



-io- 1 
icr 1 



(31a) 
(31b) 

(31c) 



The authors of Ref. [14] use the matrix 7 1 for the 
"radial" part of the Dirac equation. Specifically, near 
Eq. (21) of Ref. [14], it is stated without further ex- 
planation that a representation of the Dirac algebra is 
used where 7 1 assumes a particularly simple form, pro- 
portional to the expression given in Eq. (31a). Indeed, 
such representations exist, as we show in Eq. (31), thus 
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leading to a ramification of the somewhat ad hoc state- 
ment made in Ref. [14]. It is easy to verify that the re- 
lations {7^,7"} = 2<f !/ = 2 diag(l, -1,-1,-1) arc ful- 
filled. Here, the a = (cr 1 , ct 2 , cr 3 ) are the (2 (g> 2)-Pauli 
spin matrices, and 3.2x2 denotes the (2 ® 2) unit matrix. 



B. Eddington's Reparameterization 

In Sec. 43 of Ref. [24], Eddington has pointed out 
that a coordinate transformation exists which converts 
the Schwarzschild metric into spatially isotropic form. It 
reads as follows, 



With these coefficients, using computer algebra [25], it 
is easy to evaluate all Christoffel symbols and to establish 
that 



7 M rv 

G(r) 



7 ■ r 



w(r) G(r) . 



2v'(r) +w'(r) 
2 v(r) w(r) 



(40a) 
(40b) 



This result has been verified by us both using the repre- 
sentation (24) as well as the representation given in (25a), 
for the metric (37). 



After this transformation, the invariant line element (29) 
becomes 



C. Reduction to Radial Equation 

In our further analysis, we use the Hamiltonian form 
of the gravitationally coupled Dirac equation, 



(An + r.Y 



(33) 



1 + — ) (drj + rjdO 2 + r? sin 2 9 dip 2 
Ar\ 1 



Using this isotropic form of the metric, we can now trans- 
form the spatial part to Cartesian coordinates, 



dg2 = (An-r s ) 



(An+r s y 



+ (dxl+dyf + dzj), 

(34) 



where x\ = t\ sin# cosy, y\ = r\ sin# soup, and Z\ 
r\ cos 6. We now redefine 



n -> r, x\ ->■ x, Ux —¥ y, Zt —¥ z, 



r = \J x 2 + y 2 + z 2 . 
Furthermore, we define v(r) and w(r) as follows, 



Ar — r s ( r s 

w { r ) = z ? v \ r ) = 1 + ~r 

w Ar + r s w V Ar 



(35a) 
(35b) 

(36) 



The transformed (according to Ref. [24]) Schwarzschild 
metric can now be written in the following form, 



1 //// 



= di&g(w 2 (r),-v 2 {r),-v 2 (r) 7 -v 2 {r)) . (37) 



The considerations below are valid for a general form (37) 
of the metric and not tied to the specific form given in 
Eq. (36). The vierbein coefficients are given as 



bo !3 = b p ° = So" w(r) . 



a = a 



w(r) 



W = k 3 v(r) , (38) 
(39) 



v(r) 



where i,j = 1, 2, 3 and 8 a ^ and S a p are Kronecker sym- 
bols (i.e., equal to one if the indices are equal, otherwise 
zero). 



i(7 ) 2 dti>= + V, (41) 

where p = — i d/df with r — (x,y,z), and the expres- 
sion "Hamiltonian form" is used in analogy with flat- 
space form. Namely, in flat space, the expression on the 
left-hand side simply represents the "noncovariant time- 
evolution operator" because (7 ) — > (7 ) = 1.4x4 and 
H — > idt (see Refs. [26-29]). In curved space, with the 

metric given in Eq. (37), we have (7 ) = l 4x4 . 

In the following, we use the Dirac matrices in the Dirac 
representation, 



-0 _ / 3-2x2 

" 1 n -12X2 



~2 

7 = 



- 

a 2 
~a 2 



~5 . ~0 ~1 ~2 ~3 
7 =17777 







a'' 


t 


= (- 








1 2 X2 \ 




2x2 


) ■ 





(42b) 
(42c) 



In the form (41), the gravitationally coupled Dirac equa- 
tion allows a solution of the standard form [26-29] 



f(r) X^M 



cxp(— \Et) . 



(43) 



where the <P>) [sometimes denoted as x£(£>, p>)] are 

the standard spin-angular functions [26-29]. They have 
the property 



(44) 



We recall that the eigenvalues of the operator K = a-L+1 
are — yc [see the text following Eq. (2.9) of Ref. [27]]. In 
the text following Eq. (18) of Ref. [14], the eigenvalues is 
assumed to be +>r, an apparent typographical error. It is 
extremely instructive to write the Hamiltonian form (41), 
using the ansatz (43), in terms of (2 ® 2) spin matrices, 
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w 2 (r) 



id t tp{r) 



w(r) 



a - r ( . d . a ■ L . \ \ 

u(r) y Or r I 



3 -f I . d a- L 

-i— +i iG(r) 



l t>(r) ^ dr r 



w(r) 



ig(r) X-xn(r) 



( 



v(r) \or r ) w{r) 



v(r) \0r r J w{r) 



1/9 1 - ^ \ m „. , 

— rr 7T + + G(r) g(r) + — - / r 

w(r) \ur r / w(r) 

1 /3 h+1 „. A „, , m 
v(r) \or r J w(r) 



w 2 (r) I • i \ r~\ 



(45) 



Here, we have used the relation a ■ f x»c/i(^) = — X->t/i(^)) 
which can be found in Eq. (7.2.5.23) of Ref. [30], where 
f = r/\r\ is the position unit vector. The radial equations 
are thus given as 

o- + + G{r) g(r) = -±-L m - / r , 

Or r J w{r) \ w(r) J 

(46a) 

7T + + G(r) f(r) = [m + — - <? r . 

or r / w(r) \ w(r) J 

(46b) 

An important symmetry property of Eq. (46) is given by 
its invariance under the simultaneous replacements 



E o -E, 



f(r) o ff(r) 



(47) 



So, if is an eigenvalue of the gravitationally coupled 
Dirac equation, so is —E. Invoking the reinterpretation 
principle [31-33] and interpreting the negative energy 
— E < as +E > for antiparticles (which propagate 
"into the past" ) , we find that the spectrum of the gravi- 
tationally coupled Dirac Hamiltonian is exactly the same 
for particles and antiparticles. This important finding 
is true for any space-time metric of the form given in 
Eq. (37) and not necessarily tied to the Schwarzschild 
geometry. 

Let us now establish the connection to the flat- 
space result given in Ref. [27]. Specifically, Eqs. (2.12a) 
and Eq. (2.12b) of Ref. [27], using the identity 



r 1 d r (rh(r)) = d r h(r) + r 1 h(r), for h 
Za — > 0, can be written as 



1 and 



d 1-x 



dr r 
d 1 + h 
dr r 



g(r) = (m - E) f(r) , (48a) 
f(r)) = (m + E)g(r), (48b) 



where we have used the form (43) for the wave func- 
tion. These equations therefore become identical to our 
Eq. (46) in the limit v(r) 1, w(r) — > 1. 

The symmetry property (47) is physically tied to the 
reinterpretation principle which is very well known in 
the particle physics community [17-20] but less well 
known in the general relativity community. Some re- 
marks are therefore in order. We consider a space-time 
interval Ace = (At, Ar), and a scalar product k ■ Ax = 
\E\ At — k- Ar (with At > 0). The antiparticle amplitude 
cxp(i k ■ Ax) then is proportional to 



j\E\ At-ifc-Ar 



(-At)+ife-(-Ar) 



(49) 



where —At > and one can thus reinterpret the an- 
tiparticle trajectory, initially propagating "into the past" 
(advanced contribution to the Feynman propagator) and 
along the distance interval Ar ("from point a to point 
6"), as a positive-energy trajectory with energy +E, cov- 
ering the inverse space-time interval (—At > 0, — Ar) 
i.e., propagating into the future with four-momentum 
k), i.e., "from point b to point a." Applied to gravi- 
tational interactions, the currently available accepted in- 
terpretation based on particle physics principles [17-20] 
therefore dictates that "an antiparticle falls upward in 
the gravitational field, but backward in time, and with 
the same modulus of the kinetic energy as the corre- 
sponding particle." Therefore, after reinterpretation, the 
formalism of the gravitationally coupled Dirac equation 
predicts that antiparticles and particles receive exactly 
the same energy perturbations in a gravitational field, at 
least within space-time geometries that have the general 
form (37). This important result generalizes Eq. (1) to 
the rclativistic quantum domain. 
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IV. CONCLUSIONS 

In this article, we reexamine the gravitationally cou- 
pled Dirac equation in Sec. II, explaining a number of 
aspects of the derivation in greater detail. In particu- 
lar, we show that the condition (20) follows naturally as 
a consequence of the fundamental anticommutator prop- 
erty of the curved-space Dirac matrices (3), together with 
the known fact that the covariant derivative of the met- 
ric tensor has to vanish [10]. Under these assumptions, 
the covariant derivative of the curved-space Dirac matrix 
7^(x) also has to vanish, and the condition (20) follows as 
a consequence of the ansatz (16) for the covariant deriva- 
tive of a spinor, together with the fundamental commu- 
tator property (17). The symmetrization of the covariant 
action of the Dirac field given in Eq. (18) then becomes 
possible, in analogy to the flat-space action Eq. (4). Fur- 
thermore, under a proper definition of the local spinor 
Lorentz transformation (13), expressed in space-time co- 
ordinates, the local Dirac matrices 7 M = 7 M (x) are shape- 
invariant, as shown in Eq. (14). For a general metric 
Jf v = Tj^ix), we find the vierbein representation (24) of 
the affine connection matrices Y p = Y p (x) which differs 
from the result given previously in Eq. (8) of Ref. [12] by 
a factor j. With the additional prefactor, the result given 
in Eq. (24) then is in agreement with the result for the 
affine connection matrices given in Eq. (9) of Ref. [14]. 
In "West-Coast" conventions for the metric, the gravita- 
tionally Dirac equation reads as (i7 M V M — m) ip(x) = 
[see Eq. (19)], as opposed to the "East-Coast" form 
Of V„ + m) V>(x) = [see Refs. [12, 14]]. 

The gravitationally coupled radial Dirac equation 
given in (46) for a Schwarzschild-type metric (37) de- 
scribes the coupling of a particle (and corresponding 
antiparticle) to the gravitational field of a planet. Our 
Eq. (46), in appropriate limits, is in agreement with the 
fundamental properties of upper and lower components 
describing particles and antiparticles at rest (E — > m and 
E — > — m), if we additionally consider the limit of flat 
space-time [v(r) — > 1, w(r) — > 1]. This limit is explored 
easily, starting from Eq. (2.12) of Ref. [27]. 

The symmetry E f-> — E, f 4> g, x -f-» — x given in 
Eq. (47) implies that the quantum states of spin- 1/2 an- 
tiparticles, in the gravitational field of the Earth, have 
exactly the same exactly spectrum as those of the corre- 
sponding particles, including all relativistic corrections of 
motion. Therefore, this statement also holds for superpo- 
sitions of quantum states, including those which describe 
a wave packet evolving along a classical trajectory (these 
states would other be known as "coherent" or "Glauber" 
states). Our quantum-theory based findings go beyond 
the simple statement that particle and antiparticle tra- 
jectories in curved space-time are the same on the clas- 
sical level [see Eq. (1)]. 

Finally, let us include a remark regarding the validity 
of the gravitationally coupled Dirac equation for antipar- 
ticles. One might contemplate if antiparticles should be 
described by a different equation in the context of gravity 



than particles, but by the same equation in the context 
of electromagnetism (electromagnetically coupled Dirac 
equation). In this case, the gravitationally coupled Dirac 
equation (19) would only describe particles, not antipar- 
ticles, even if it admits negative-energy solutions. How- 
ever, in this case one gets into trouble in the limit of a 
vanishing gravitational interaction, in which case space 
becomes flat. This is because the Dirac equation is known 
to describe antiparticles very well in this limit [17-20]. At 
least, this concept is used in all perturbative QED cal- 
culations, including the notoriously difficult bound-state 
problems [34]. If we conjecture that the flat-space limit 
is smooth, then the gravitationally coupled Dirac equa- 
tion (19) must remain valid for both particles and an- 
tiparticles. Brill and Wheeler [12], Boulwarc [13], and Sof- 
fcl, Miiller, and Grciner [14] all used the same methods for 
deriving the coupled Dirac equation (we here attempt to 
resolve some discrepancies found in the literature regard- 
ing the final steps in the derivation) . The gravitationally 
coupled Dirac equation (19) involves (4(g>4) matrices and 
allows for two positive-energy solutions, which are natu- 
rally interpreted as particles, and two negative-energy so- 
lutions, which are naturally interpreted as antiparticles, 
according to usual practice in particle physics [17-20]. We 
may thus assume that the gravitationally coupled Dirac 
equation given in Eq. (19) should be valid for particles 
and antiparticles, simultaneously. We know for sure that 
the corresponding variant of the equation in flat space 
describes particles and antiparticles simultaneously, as 
described in detail, e.g., in Chap. 2 of Ref. [19]. 

Despite our theoretical considerations, it is still of ut- 
most value to the scientific community to carry out the 
planned antimatter gravity experiments [3-8]. We con- 
clude with two remarks: (i) The "inertial mass term" 
in the sense of the equivalence principle enters the 
Schrodingcr equation for free and bound electrons, e.g., 
for a bound electron in a hydrogen atom. According to ex- 
perimental evidence, inertial and gravitational mass are 
the the same for atoms (such as atomic hydrogen), which 
is composed of spin-1/2 particles (electrons and protons), 
Therefore, the action principle 6S = [see Eq. (11)] 
provides for a solid basis of the discussion of relativis- 
tic quantum effects in gravitational coupling, with the 
m term entering the equation being equal to the iner- 
tial (gravitational) mass, (ii) Our investigations suggest 
that any conceivable differences of the gravitational cou- 
pling of particles and antiparticles should be assigned to 
a "fifth force," not to any conceivable "modifications of 
the gravitational mass" of antiparticles versus particles. 
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